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INVESTIGATION OF SUPERSONIC TURBULEITT BOUNDARY- LAYEl^ SEPARATION 
ON A COMPRESSION RAMP BY Al'I INTEGRAL METHOD 

D. K. Patel and K. R. Czarnecki 


SU!-HMARY 


An investigation vas made to determine the feasibility of using a 
boiindary-layer integral method to stud^r the separation of a turbulent 
bo-undary layer on a two-dimensional ramp at supersonic speeds. The .numeri- 
cal calculations were made for a free-stream Mach number of 3, a Reynolds 
number of 10 million, and over a ramp-angle range from 0° to 30°. 

For ramp angles where no flow separation was indicated, the theoretical 
calculations were in reasonable agreement with experimental data except 
for a somewhat belated rise in the pressure. For the larger ramp angles, 
where separation was present, the investigation produced results that 
were not in agreement with experiment or with results calculated by time- 
dependent Navier-Stokes methods apparently because no provision had been 
made for a proper shock-boiondary layer interaction where strong normal 
pressure gradients are induced within the boiindary layer under the shock 
independent of surface curvatiu'e effects. Within the limits of the 
calculations, the effects of changes in the size of the fillet, the 
omission of curvature terras, and changes in eddy- viscosity model were neg- 
ligible. 


INTRODUCTION 

With the development of large high-speed computers , one of the aero- 
dynamic problems currently receiving considerable theoretical attention 
is that of the separation of a turbulent boundai^'" layer in a two- 
dimensional compression corner at supersonic and hypersonic speeds. 



Although progress is being made in theoretically predicting separation 
characteristics which agree reasonably well with experiment (see refs. 1 
and 2) many aspects of the problem remain unsolved. In particular, for 
combinations of free-streara Mach number and ramp angle for which there is 
no shock detachment in inviscid flow and if the boundary-layer thickness 
is small relative to ramp chord the length of the separation region is 
dependent upon boundaiy characteristics alone. For combinations of Mach 
nijmber and ramp angle which result in the shock detachm.ent angles being 
exceeded the separation length is dependent primarily upon the maximxmi 
height of the ramp if the boundary layer is relatively thin and on both 
the ramp height and boundary-laj^er characteristics if the boundary layer is 
relatively thick, in close analogy to the for'jard-facing step flow 
separation (ref. 3). The transition from one regime to the other needs 
further investigation. This investigation can best be done by a param- 
etric study. Unfortmately , the time -depen dent IJavier-Stokes methods 
which appear to give reasonable results (refs. 1 and 2) are rather 
expensive and lengthy timewise for such an investigation. Consequently , 
it appeared desirable to determine if it was feasible to make such a studj'' 
using a boundary- layer integral method where the inherent simplifying 
assumptions generally eliminate significant am.oiuits of computer time. 

The boundary-layer integral method used in this investigation is 

based on the method devised by Kuhn and Nielson (refs, h and 5) and 

extended by Tai to include the strong-interaction case (ref. 6). The 

same velocity profiles and eddy- viscosity models utilized by the afore- 
/ 

mentioned investigators were used, but the basic boundary-layer equations 
were modified to include surface curvature terms so that a singularity 
at the intersection of the ramp with the flat plate could be eliminated 
by incorporating a fillet. External- flow characteristics were determined 
assuming a Prandtl-Meyer supersonic flow. No investigation was made of 
the case where the external flow behind the boundary- layer reattachment 
point was subsonic or of the characteristics of separated-flow regions 
because of a major deficiency in the general method. Both the weak- 
and strong-interaction methods of Tai were investigated along with 
many modifications. 
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Nijmerical calculations vere made for a free-stream Mach number of 
3, over a ramp angle range from 0° to 30°, at a Reynolds nmber, based 
on the flat-plate length to the ramp leading edge, of 10 million. The 
fillet between the ramp and flat plate was varied from 1.5 to 6 percent 
of the flat-plate length to the ramp and its shape was varied. Most of 
the solutions were obtained by a point-iteration downstream-marching 
procedure but some global iteration solutions were also obtained for 
comparison. 


KOfiENCLATURE 


a 


c. 

1 

D 


f(n) 

F(n) 


k 


I 


m 

M 


+ 


n 


speed of sound 

local skin- friction coefficient. 


A 

T 

w 


, A A 

l/2p u 


coefficients of polynomial, equation A-1 

distance from plate-ramp comer to beginning or end of 
fillet , figure 10 

functions defined in equations C-1 through C-9 
weighting function, equation 22 
function defined in equation C-13 
constant in equation 35 

distance from leading edge of flat plate to plate-ramp 
comer 

constant in equations 3^ and 35 
Mach number 

dimensionless normal distance from plate-ramp surface , ^ 

V 
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p p P-> 


exponential constant in equation 3 


dimensionless static pressure ?- 

P 




functions defined in equations C-10 througii C-12 


dimensionless radius of curvature , ^ 


Reynolds nximber , 


A / 

p u 

'^OO CO 


s ,n 


dimensionless orthogonal curvilinear coordinates measured along 

s n 


and normal, to plate-ramp surface ; ^ £ 

time 


dimensionless temperature. 


A 

T 




dimensionless velocity components aiong and normal to plate- 

A A 

ramp surface ; 


A ’ A 

u u 
00 00 


dimensionless wake velocity. 


A 

A 

U 


A 

U 


dimensionless friction velocity, - — 

u 


friction velocity. 



V 


dimensionless resultant velocity 


f~2 2 

, Wu + V 
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x,y 


cartesian coordinates 
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a 

6 


^ 1 ’^ 2 ’^3 


y 

6 


A 

e 


ramp angle 

eddy viscosity, defined in equations 27-30 
function defined in equation 3 
eddy-viscosity parameters, equation 19 

specific-heat ratio 

A 

dimensionless bomdary-layer total thickness, — 

^ A* 

dimensionless boundary-layer displacement thickness — 

I 

eddy viscosity 

A 

0 

dimensionless bo\mdary-layer momentum thickness 

SL 

boundary-layer flow angle defined by equation 38 or 
equation Ul 

exponent denoting temperature-viscosity relationship, 
equation 3 

viscosity 


V 


Prandtl-Meyer angle 


A 

V 


kinematic viscosity 


density 


exponent defined in equation 22 

A 

T 


dimensionless shear stress. 


A /\ 2 
P u 

^OO CO 


Subscripts 


e 


compressible flow 
edge of boundary layer 
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m median 

t total or stagnation 

w wall 

°° free-stream 

Superscripts 

' turb\ilent fluctuating velocities 

time averaged mean quantities 
dimensional quantities 

Note: All quantities without the subscript c are in the incompressible- 

flow plane. 


THEORETICAL MODEL 

Details of the theoretical model used in this investigation are 
shown in figure 1. A two-dimensional supersonic turbulent boundary layer 
on a flat plate approaches and then is deflected onto the ramp surface at 
corner C. If the ramp angle is sufficiently small the boundary layer 
remains attached during the turn and ens\iing encounter with a strong 
adverse pressure gradient. In order to eliminate a singularity in the 
computations a fillet is incorporated into the comer, the size and 
shape of which can be varied to influence the shape of the pressure rise 
if necessary. (For details of fillet design see Appendix A.) Thus, the 
usual thin comer shock is essentially replaced by a set of strong coalescing 
compression waves and the expected abrupt pressure jump by a steep but 
finite adverse pressure gradient. Inasmuch as boundary-layer separation 
might be expected to first occur on the fillet as the ramp angle is 
sufficiently increased this approach requires that the basic boundary- 
layer equations retain curvature terms to r~0(6). The inviscid flow 
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above the boundary layer was assumed to be governed by the Prandtl-Meyer 
relations. Calculations indicated that the maximum expected entropy 
changes should be small and their influence on the flow characteristics 
negligible for the range of conditions that were to be explored in this 
investigation. 

For greater simplicity the viscous boimdary-layer calculations were 
made in the constant property or incompressible-flow plane with the 
transition from or to the compressible boiindary layer being made via 
Stewartson's transformation (ref. ?)• The interaction between the 
boundary layer and inviscid exterior flow w^as investigated in the 
compressible-flow plane by two methods. In the weak interaction 
(see ref. 6) the angle of the inviscid stream was matched to the angle 
generated by the growth of the displacement thickness of the boundary 
layer and the velocity at the edge of the bo\mdary layer was then matched 
to that of the stream. In the strong interaction the velocity and flow 
angle of the inviscid flow were matched to the velocity and streamline 
angle at the edge of the boundary layer. 


INVISCID FLOW 

The pressure, density, Mach number, etc., characteristics of the 
inviscid flow over the boundary layer were determined by Prandtl-Meyer 
flow. Preliminary calculations based on experimental indications 
revealed that the entropy changes that might occur in actual physicsil 
flows were relatively small for the Mach number range of the investigation 
and that their omission in the present theoretical approach would have 
no significant effect on the general flow characteristics to be calculated. 
Actually, the cadculations were made by using the Prandtl-Meyer flow 
tables of WACA 1135 (ref. 8) with the Prandtl-Meyer angle V being 
used as the entering variable to make the interpolations. The angle V 
is foiand from either- the slope of the boundary-layer displacement thick- 
ness or the angle of the flow at the outer edge of the boundary layer 
in the compressible- flow plane. 
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VISCOUS FLOW 


Stewartson Trans fonnat ion 

As mentioned earlier the hoimdary-layer calciilations were made in 
the incompressible- flow plane whereas the interaction between the 
inviscid and viscous flows were carried out in the compressible-flow 
plane. Ibe transformation from one plane to the other was accomplished 
by means of the Stewartson transformation (ref. T)- 
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Additional transformation eq^uations can be derived directly from the 
definitions of the pertinent quantities and the Stewartson transformation 
or from the definitions, transformation and Crocco's relationship between 
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locaJL temperatures and velocities within a boundary layer having a constant 
static pressure and total temperature across the layer. Some of the more 
useful relationships are 



The transformation equation derived from the Stewartson transformation 
for the local, skin- friction coefficients was not considered reliable 
because it depended upon the temperatures or densities at the wall 
whereas the T' methods indicate a mean boundary-layer temperature or 
density approach is more reliable. (See ref. 9 , for example.) For this 
investigation therefore the assumption was made for simplicity that such an 
equation derived in reference 10 would be adequate 
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(3) 


= B (1 - (1 + X) H) 


where 


* 

B 


(1 + 


x.-t 1 


) 

e ,c 


-1 


X =0.76 


N 



Governing Equations 

The basic boimdary-layer equations were derived from the two-dimensional 
incompressible Navier-Stokes equations in the curvilinear orthogonal coordi- 
nate system (ref. 11, page 98): 
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( 6 ) 


In the analysis the steady-state case is considered, thus, the — — terms 

3 “t 

hecome zero in equations (5) and (6). 

To convert the above equations to dimensionless form the following 
dimensionless ratios were adopted 
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A 
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A A ^ 

Here P = P<jo incompressible flow. % is the distance between the 

flat-plate leading edge and the ramp comer. 

Thus, equations (U) through (6) become: 


r a u a V V 

r + n3s an r + n 


0 


( 7 ) 
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r 3 u 9 u V u 

U — + V -T — + 


r+n 9s 9n r+n 


r _Le. 
r + n 9 s 


R 


9 9 ^ 1 9 u 
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I . \2 ^2 -2 r + n9n / .\2 

(r + n) 9s 9n (r + n) 


2r 9 V 


r 9 r r n 9 r 9 u 

, v2 9s , >3 9s^^, . v3 9s 9s 

(r + n) (r + n) (r + n) 


( 8 ) 


r 9 V 9 V _ _u 

r + n 9 s 9 n 


= _ JLi+ L. 


r + n 


9 n R 


9 V 


9n 


2 r ^ + 1 9 V ^ r. 


9 V 


/.>2 9s r + n9n r . -s ^ ^ 

(r + n) (r + n) 9s (r + n) 


r ^ r n ^ 9 v 

/ n3 9s'^ / ^3 9s 9s 

(r + n ) (r + n ) 


(9) 


Equations (7) through (9) are reduced to houndary-layer equations by 

the standard order of magnitude analysis iised for shockless boundary layers. 

9 

This procedure does not take into account that under the shock the -r — 

g 9 s 

terms are on the same order of magnitude as the — terms and, consequently, 

d n 

large normal pressure gradients may be induced by terms other than those 

9 V 

involving centrifugal forces. In particiolar, the '^“ 5 ^ term is probably 
very important in generating normal pressure gradients and having a strong 
influence on the boundary-layer separation characteristics under a shock. 
Still, it is of interest to determine whether the significant characteristics 
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of the separation can be reproduced by inducing separation by means of the 
adverse gradients generated by a fillet and making the calculations by means 
of a second-order boundary-layer theory. Experimental results indicate that 
the pressure rises associated with shock- induced separation occur in one or 
two boundary-layer thicknesses (see ref. 12, for example); thus, the fillet 
must be relatively small. Estimates indicate that curvature terms of 
r = 0(6) must be retained in the equations. The results are: 


Continuity: 


r 9 u ^ 9 V ^ V 

r + n9s 9n r + n 
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(10) 


s -Momentum: 
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r + n'^9s "'^9n r + n r + n9s 
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« 9 n 
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( 11 ) 


n-Momentum: 


9 p _ u 
9 n r + n 


(12) 


To convert the equations (lO) through (12) to turbulent-flow equations, 

let 
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u = u + u' 

V = V + v' 
p = p + p' 

and after substitution of the above vaJ-ues in equations (lO) through (l2) 
and after the equations are mean time averaged, the following equations are 
obtained 




r+n3s 9n r+n 
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9n r + n r + n 


(15) 


In the shock interaction region the magnitude and characteristics of the 

term u'^ is not well known (ref. 13). Thus, for this analysis the term 
2 

u' is not taken into account. 

In order to eliminate the Reynolds shear stress, an eddy viscosity is 
introduced as 
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Thus, equations {ih) and (15) become 


s -Momentum: 


r - 3 u - 3 u 
— 7 — u -r— + V -r — 
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and 


n- Momentum: 
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(18) 


vhere 


B. = 


1 + 1 


2e 




^2 = ""i-="^i-" r 


83 = 1 
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(19) 


Equations (13), (IT), and (I 8 ) are the governing equations for the viscous 
flow. The boundary conditions are u = o, v = v^ at the surface and u = ^^(s) 
at the edge of the boundary layer for the weaic interaction and u = o , v = v^ 
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V 2 2 

at the edge of the boundary layer for 

the strong interaction, represents a surface injection or suction velocity 

if one is desired. 


Integral Method 

The desired integral equation is obtained by combining the continuity, 
s- and n-momentiam equations (eqs. (l3, ( 17 )» and. (l8)) in a manner similar 
to that of Tai (ref. 6) and Ku hn and Nielsen (refs. 1+ and 5)* Equation (13) 
is solved for the normal velocity 



( 20 ) 


and the n-momentum equation (eq. (l8)) is integrated vith respect to n 

to remove the n dependency and differentiated with respect to s to 

3p 

obtain the pressure gradient as 
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(21) 


where n is a dummy variable of integration. Equations (20) and (2l) 
are substituted into equation (IT) and the resulting equation is multiplied 
by a weighting function 


f(n) = n^; a = 0, 1 


(22) 


and integrated completely across the boundary layer to yield the desired 
integral equation as 


16 


ORIGINAL PAGE IS 
OS EOOR QUALITY 



OMGINAC page is 
Oi^ POOR quality 


6/_ n__ _ n_ 

f ) “ 3 u f u 3 u , / 3 u . u \ r 3 u 

Jj^3s ^Jr + n3s^’^ \3n r + n/j3s‘^’^ 

r\ * A A 


^ r + n 3^ 

3 s r R 


— ( 3 . 4^) + ^ ~ 


3n 1 3n r + n 3n , ^2 

(r + n) 


r -2 

-I* 


n) 


3 r 
2 3s 
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In equations (22) and (23), CT = 0, 1 corresponds to the momentum and moment- 
of-moment\jm equations, respectively. 

The two equations derived from equation (23) using 0=0 and 1 are 
sufficient to obtain solutions for the case of the weak interaction. For 
the case of the strong interaction, which might be expected to apply in 
the neighborhood of the fillet between the flat plate and ramp, an additional 
equation is required to obtain a solution. This equation is derived by 
integrating the continuity equation (13) completely across the boundary 
layer (essentially extending the limit n in equation (20) to 6) to get 



r + 6 - 

V + 

r e 


V 
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(2U) 


Velocity Profile 

The integration of equations (23) and (24) is accomplished by 
assuming a velocity profile that has the potential for indicating separated 
(reversed flow) characteristics in strong adverse pressure gradients. A 
profile that had satisfactory characteristics in subsonic and transonic 
flow on a convex surface was developed by Kuhn and Nielsen (ref. 4) and is 
utilized in this analysis 
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= u [2.5 5,n (1 + n"^) + 5-1 - (3.387 n'*’ + 5.l)e ] 
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/ n \ 

+ 0.5 Ug 

1 - cos (tt j ] 



(25) 


where 




is the wake velocity and 


u 

T 


is the friction velocity 


u 

T 




(26) 


In eciuation (25), the first term on the right-hand side represents the 
inner part of the velocity profile , consisting of a laminar sublayer 
whose function is to allow the no slip boundary conditions to be met at 
the wall and of the law-of-the-wall function. The second tenn represents 
the wake portion or outer part of the boundary layer by a cosine function 
that is easier to handle in the calculations than the more accurate error 
function representation. 


Eddy Viscosity 

The eddy viscosity model used in the present analysis is identical 
to that of Tai (ref. 6) which is similar to that of Kuhn and Melsen 
(ref. 5). 

The expressions for the eddy viscosity are given as follows: 

For attached flow, inner layer 


8=1+ 0.0533 
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1 + o.Ui 
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(27) 


For attached flow, outer layer 
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0.013 + 0.0038 e 




dp 1 
ds 15 
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6 R 


03 


(28) 


For separated flow, inner layer 


6=1+ 0.018 u n R 
e e 




(29) 


For separated flow, outer layer 


0.013 u 6 R 


6 = 


1 + 5 


•<f) 


(30) 


System of Equations for Strong Interaction 


Equation (25) is substituted into equations (23) and (24) with u^ 

p 

being eliminated by evaluating u at n = 6. Details as to how this is 

actually accomplished are given in Appendix B. -Three ordinary differential 

T d£ *^^w 

equations are finally obtained with \mknowns , and -r — irhich can 

ds ds ds 

be written in the matrix form 
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functions 
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du 
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Q. 


Q, 


(31) 


and Q. are presented in Appendix C. TJote that 
1 J J 

because of the existence of surface curvature effects the substitution of 
equation (21) into equation (l?) makes it more convenient to make the wall 

^ the third unknown in the present investigation^ - 

S_J L _J_ __1 3 T 

3 s 


pressure gradient 
rather than the velocity gradient at the edge of the boundary layer 


as was done by Tai who had no such terms. Furthermore, the experimental 
pressure gradients with which theoiy is compared are wall pressure gradients. 

Systems of Equations for Weak Interactions 


Most of the calculations for the weak interaction case were made by. 
an approach similar to that of Tai 
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du n 
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where E. . 

ij 


and Q . 
J 


are identical to those in equation (3l). As u 


tended toward zero near separation the coefficient determinant also tended 
toward zero with the result that the calculations diverged. For some of 
these conditions the solution algoritiim was chanp.;ed to 
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( 33 ) 


where an initial value of u^ was prescribed by a parabolic arc extrapo- 

T 

lation from previous stations so that the derivative ^ could be an 

ds 

input. Tlie solutions of equation (33) were then iterated with the 
Prandtl-Heyer inviscid flow, with continuous updatix'.g of u^ , until the 
velocities calculated for the outer edge of the boundaiy layer by the 
viscous equations (eq. (33)) matched the velocities computed for the 
Prandtl- Meyer flow. 


Initial Conditions 

Initial values for the boundary-layer variables were derived from 
Schlichting’ s equations for an incompressible flat plate boundary layer 
(ref.' 11, pp. ^33 and 468) modified for taking account of a pressure 
gradient 


and 
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where 


m = I4, k = 0.0128 

These equations were derived assuming a l/7th power velocity profile , hence , 

are not compatihle with the velocity profile used in this investigation 

- « 

(eq. 25). Consequently, u_^ and 6 from equations (3^) and (35) were 
iterated with the velocity profile of equation (25) to obtain new values of 
boundar;'’- layer variables that were compatible both with that velocity profile 
and the pressure gradient existing at any particular iteration. 

COUPLING OF INVISCID MD VISCOUS FLOV/S 

In the case of the strong interaction the requirements for the coupling 

of the inviscid and viscous flows are that both the magnitudes and directions 

of the velocities of the two flows match one another at the outer edge of 

the boundary layer. Because of entrainment of air from the inviscid into 

the viscous flow as the bomdary layer grows along the surface , the slopes 

of the streamlines at the edge of the boundary layer are smaller than the 

d.5 

slopes of the surface of the boundary layer — . The relation between the 
slopes is readily found by the use of the continuity equation to be 


d6 . 
m— = tan 
as 


G + 


P u. 
e e 


ds 


[p u (6-6 )] 
e e 


which can be written as 


(3b) 


= tan G + (6 - 6 ) -^ in (p u ) 
ds ds e e 


( 37 ) 
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and solved for the streamline angle 


£ = tan ^ (6-6 ) £n (p u )] 

as as e e 


( 38 ) 


The normal velocity in the inviscid flow at the edge of the boundary layer 
is found from 


V . = V. sin 0 

e,inv inv 


( 39 ) 


and the edge velocity in the equation defining the assumed velocity 

distribution within the boundary layer, equation (25), is thus related to 
the inviscid flow by 


u 

e 


/ 2 
(v, 
inv 


e ,inv 


iho) 


For the case of the weaii interaction the slopes of the inviscid stream-^ 

d(S 

lines are matched to the growth in boundary-layer displacement thickness 
and equation (38) reduces to 


O + -1 V 


( 41 ) 


Also, the assumption is usually made that 


V . is negligible, hence 
e ,inv 
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U = V. 

e inv 


( 1 * 2 ) 


Because of divergence problems that were encoimtered in the strong inter- 
action calculations some weak interaction solutions were obtained both 
with and without this latter assumption. 

COMPUTATIONAL METHODS 

The numerical integrations of the viscous-flow equations in the s 
direction were carried out by means of a standard fourth-order Runge-Kutta 
method for both the point-iteration and global-iteration techniques. In 
the point-iteration technique the interaction between the viscous and inviscid 
flows was iterated at a constant s station until the solution converged, 
broke down or a prescribed niimber of iterations were completed. In the 
global-iteration method the boundsiry-layer calc\ilations were executed over 
the complete s domain before reevaLuation of the new inviscid pressure 
gradients and recalculation of the boundary layer. The computations were 
usually begun at s = 0.20 and continued downstream until the calculations 
broke down or an s of 1.5 was reached. Relaxation factors of varying 
magnitude were used on one or more of the dependent variables and various 
types of smoothing procedures were incorporated into the global- iteration 
calculations. All calculations presented herein were made with equal 

to zero. For most calculations the Stewartson s transformation was ommitted 
because its effects were found to be generally negligible and its presence 
increased the complexity of the calculations. 

As the calciilations proceeded toward a separation point in the weak- 

interaction point-iteration method the solutions obtained by Tai’s method, 
du 

where and — were the imknown dependent variables , began to diverge 

ds ds 

because the determinant for the equations tended toward zero. For some cases 
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the equations vere revised so that 

became the \mknowns . Ttie gradient 
extrapolating parabolically from the previous stations and calculating 

a three-point backward derivative using the extrapolated u^ as the third 
point. During the iterations u^ was continuously modified to require 
that the boundary-layer edge velocities in the inviscid and viscous flow 
were identical within prescribed limits. In many ways this type of calcula- 
tion scheme seemed to be more naturally suited to supersonic flow problems 
than the approach used by Tai (ref. 6). 


du dp 

T , do , -^w 

- — was prescribed and and — 

ds ds ds 

du 

X 

- — was initially prescribed by 


RESULTS MD DISCUSSION 

General Comments on Stability and 
Convergence of Calculations 

Before any final results of the calculations are presented it appears 
desirable to point out a few interesting or unusual characteristics that 
were encountered in making the computations that are not apparent from the 
final plots themselves. These are: the extreme instability of the calcula- 

tions ; the unusual values of relaxation factors that were required to allow 
the calculations to proceed; the ineffectiveness of the relaxation coefficients 
in improving convergence; and, finally, the inability to obtain any useful 
results using the strong-interaction or global- interact ion techniques. (See 
ref. lU for similar problems in instability in calculations in large adverse 
pressure gradients using a similar integral approach. ) 

First, there never was any problem in calculating the boundary-layer 

characteristics on the first iteration where the pressure gradient was 

assumed to be zero. As the boundary layer was iterated with the inviscid 

flow in the second and higher iterations many of the terms involved in 

the calculations tended to oscillate in value with the magnitude of the 

* 

oscillations usually being catastrophically large. As an example, 6 
would often increase so much between adjacent stations as to make the 
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so large that the inviscid-flow calculations would indicate 


. d6 

slope 

ds 

the existence of a detached shock with a subsonic inviscid flow above the 
boundary layer. This would abort the calculations because no provision 
had been made for the existence of subsonic flow this early in the investiga- 
tion. Again, calculations would often abort because 6 would become 
negative d\xring the oscillations and no provision had been made in the 
equations for the interpretation of negative values of boundary-layer total 
thickness. In general, only a few of these types of problems could be 
solved without recourse to the introduction of relaxation factors. 

The role of the relaxation factors was to introduce the interaction 
effects into the computations slowly and hence decrease the magnitude of 
the oscillations. Because of the magnitudes of the oscillations that 
existed without relaxation, nmerical values of relaxation factors an order 
of magnitude smaller than \isual had to be used. For example , in steep 
adverse pressure gradient regions but still some distance away from separa- 
tion points, relaxation factors of 0.05 to 0.2 were required in the 
weak-interaction point-iteration method on two or more of the most critical 
terms to preserve order in the calculations. These small values of relaxa- 
tion factor required about 200 iterations to permit satisfactory convergence. 
In the global-iteration scheme single relaxation factors on the order of 
0.2 to 0.5 were sufficient, but this was because smoothing procedures were 
simultaneously employed. 

A short study was made to determine the effects of relaxation factors 
on convergence. In the study no evidence was ever fo\xnd that the use of 
relaxation factors converted a divergent solution into a convergent one. 

In fact , all evidence indicated that the use of relaxation factors in 
divergent solutions merely delayed the time required to reach the point 
where the computations broke down (i.e. , larger mmiber of iterations). 

Attempts to obtain converged solutions by the strong-interaction 
technique were universally unsuccessful. In the point-iteration method 
these attempts included starting the calculations at s = 0.20 by the 
use of the usual initial procedures, starting the calculations just ahead 
of the fillet from the boundary layer determined by the weak-interaction 
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technique, and using a range of relaxation factors. Sven the use of the 
glohal-interaction method resulted in divergent solutions for the flat 
plate . 

In veak interactions converged solutions could he obtained hy the 
use of the global-iteration technique combined with a cubic-spline smoothing 
procedure. Without smoothing, oscillations quickly began to appear in areas 
such as the jimctures of the fillet with the approach surface and the ramp 
as the calculations proceeded downstream. The magnitude of these oscillations 
increased rapidly with increases in the number of iterations but without 
changes in location. Unfortimately , the smoothing procedure first controlled 
then decreased the magnitude of the oscillations by merely flattening the 
waves in both the upstream and downstrearii directions. Inasmuch as the positive 
and negative oscillations were completely out of phase with one another the 
oscillations were not properly cancelled and the converged solutions were 
obviously physically incorrect. In essence, the smoothing procedure allowed 
the errors from the critical points to seep in until the flow changes at the 
critical points were overwhelmed by the smoothing corrections, whence conver- 
gence could be achieved. 

An examination was also made of the effects of grid size on stability 
of calculations and convergence. The results indicated no effect except to 
pinpoint more sharply the station at which divergence actually began. 

Finally, the u^ initially prescribed method appears to alleviate many 
of the previously mentioned problems, but its potential was not fully 
explored. 


Effects of Various Terms or Parameters 

Some final results computed by the weak-interaction method with point- 
iteration technique are presented in figure 2. In the figure, the non- 
dimensional wall-pressure ratio p^ and the local skin- friction coefficient 
c^ are plotted as fimctions of surface distance s for a range of wedge 
angles from 0° to 30°. The fillet is based on a 9th power polynomieil and 
begins at s = 0.98 and ends at s = 1.02. Theoretical computations are 
represented by the lines ; the symbols have been included at intervals only 
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to aid in identifying the curves. Only converged data are plotted. Data 
for c^> 0.0008 were computed by the Tai method and data for c^ < O.OOO 8 
for a = 30° were computed by the u^ initially-prescribed method. Similar 
calculations for c^ < 0.0008 were not made for the other wedge angles for 
reasons that will be discussed in the next section. The ticks on the pressure 
plot indicate the inviscid-flow pressure rises expected for the various wedge 
angles using shock equations. 

The curves of figure 2 indicate that as the wedge angle increases the 
adverse pressiire gradient also increases and the local skin-friction 
coefficient decreases. Boiindary-layer separation was calculated for 
a = 30° at s =1.007 and appears to be a probability for a = 25°, 20°, 
and 15 °. llote that the pressure rises lag changes in surface curvature. 
Approximately one-half of theoretical pressure rise at a = 10° occurs 
on the flat wedge surface downstream of the fillet. A comparison of the 
a = 10° calcvilations where no separation is indicated with some experim.ental 
data obtained at a Mach number of 2.75 (ref. 15) indicates (fig. 3) that, 
although the theoretical pressure rise occiirs somewhat belatedly, it is in 
reasonable agreement with the experimental results. 

Some nondimensional velocity profiles corresponding to the wedge 
angle of 30° are shora in figiire U. They corrobrate the appearance of 
separation. Note that practically all the significant boundary-layer 
profile changes occur in an extremely short region that corresponds to 
the steepest part of the adverse pressure gradient in figure 2 . 

At this point the question arises as to the possible effects the 
fillet could have on separation. Figure 5 shows that even for the larger 
wedge angles the effects of changing the size and, hence, radii of curva- 
ture of the fillet had very little effect on either the pressure or skin- 
friction distributions. Figure 6 shows that the surface curvature terms 
contributed little or nothing despite the relative smallness of the radii 
of curvature (lowest r /6 was about 2 . 9 )- 

Finally, figure 7 indicates that sizable variations in the edd^’'- 
viscosity model were of no significance. 
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Calc\ilations With Experimental Input 


Although the results disc\issed in the previous section would appear 
to indicate the possibility of deriving a reasonable criterion for boundary- 
layer separation a closer examination reveals a serious problem. The 
data of figure 2, for example, indicate that separation occurs at wall 
pressure ratios approaching 3 or more whereas experiments such as those 
in references 12 and l6 and theoretical calculations such as those in 
reference 17 indicate that separation at = 3 occurs for pressure ratios 
less than 2. In order to test this discrepancy some calcxilations were 
made using Law's (ref. 12) experimental pressure distribution for a = 25°, 
where separation was well developed, to represent the final converged 
pressure distribution. The resul.ts are presented as figure 8. Included 
in figure 8 are the theoretical skin-friction results of Shang and Hankey 
(ref. 17 ) whose pressure distributions calculated by time-dependent 
Navier-Stokes methods are in reasonable agreement with Law's experimental 
data. 

The data of figure 8 and the discussion in reference 12 indicate that 

boundary- layer separation occurred experimentally at an s of approximately 

0.96 and a p of approximately 1.9* The calculations of Shang and Hankey 

(ref. 17 ) indicate separation at s approximately 0.95 and. p about 1.03. 

w 

Present calculations indicate separation at s = 1.07 and p^ = 4.3. 
Furthermore, the experimental separation point occurs in the first pressiire 
rise leading to the pressure plateau region whereas in the present investiga- 
tion separation occurs in the second pressure rise at a much higher pressure 
in a region where a reattachment shock might be expected. This discrepancy 
is ascribed to the failure to provide in the present investigation for 
terms which can induce normal pressure gradients withou the need of surface 
curvature. This omission becomes even more serious for the cases where the 
ramp angle is large enough to force separation to occur well ahead of the 
fillet where there is no surface curvature and hence no mechanism for 
boundary-layer separation in the present approach. As illustrated in 
figure 9, due to the strong coelescence of the compression waves into a shock 
as one proceeds away from the surface, the pressure gradients at the wall 
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are considerably smaller than those at the outer edge of the boundary layer 
or the mean gradient that is essentially utilized in the integral-equation 
calcijlations . Omission of the tenris that generate the normal pressure gradients 
also means ommission of the coelescence effects at the edge of the boundary 
layer and consequently much smaller mean axial pressure gradients than those 
that actually should exist. A similar serious problem was encomtered by 
Werle and Eertke (ref. l8) in an investigation of the supersonic wedge problem 
wherein the boundary-layer equations were solved by a finite-difference 
approach. This concept was further substantiated in the present investigation 
by making some calculations wherein the pressure gradient within the first 
pressure rise was arbitrarily steepened without modification of the plateau 
pjressure. Eoundary-layer separation was readily achieved at the desired 
location. Additional evidence that normal pressure gradients in the- boundary 
la2rer must be accounted for in integral methods is presented in reference 19 
’.fnere the boundary layer is subsonic but is thick and has significant stream- 
line curvature. 

iiodifi cat ions can be made to the integral-equation metiiod to provide 
for a more realistic shock-bovindary-layer interaction but at the present 
time the modifications will have to be either empirical or built up by 
interative inclusion of one or two more terms. From the standpoint of the 
desired end result the normal pressure gradients should not be prescribed 
a priori but should be an integral part of the interaction solution. In 
this respect, both the inviscid and viscoias flow solutions in the present 
investigation are essentially one- dimensional calcialations . It is possible 
that a practical solution of the shock-bomdary-layer problem will reqiiire 
calculations in two dimensions. In view of the major deficiency just dis- 
cussed in the present approach it was felt that further calculations without 
major modifications were not justified. 

The strong normal pressure-gradient effects in shock-boundary-layer 
interactions are apparently peculiar to turbulent boimdary layers. In 
laminar flows the pressure rises and pressure gradients, even at separation 
points, are an order of magnitude smaller and pose no major problems in 
the use of integral methods (ref. 20). 
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SUI-S’IAEY OF PtEGULTS 

An investigation was made to determine the feasibility of using a 
boundary- layer integral method to study the separation of a turbulent 
boundary layer on a two-dimensional ramp at supersonic speeds. The 
numerical calculations were made for a free-stream Mach number of 3, a 
Reynolds number of 10 million, and over a ramp angle range from 0° to 30°. 

For ramp angles where no flow separation was indicated, the theoretical 
calculations were in reasonable agreement with experimental data 
except for a somewhat belated rise in the pressure. For the larger ramp 
angles, where separation was present, the investigation produced results 
that were not in agreement with experiment or with results calculated by 
time-dependent Navier-Stokes methods apparently because no provision had 
been made for a proper shock-boundary layer interaction where strong normal 
pressure gradients are induced within the boundary layer imder the shock 
independent of surface curvature. Within the limits of the calculations 
the effects of changes in the size of the fillet, the omission of curvature 
terms, and changes in eddy-viscosity model were negligible. 
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APPENDIX A 


To avoid a singularity in the computations at the sharp corner C 
a fillet is introduced as indicated in figure 10. For best results it is 
desirable that the entry to the fillet at A and exit from the fillet at 
B be as smooth as possible and the smallest radii of ciorvature be con- 
centrated close to the comer C. These conditions can most easily be 
attained by the \ase of a high-pover polynomial wherein several of the 
higher order derivatives at A and B are set to zero. Preliminary calcula- 
tions indicated that a 9th power polynomial woxild be optimum for the calcula- 
tions. Some calc\ilations were also made with a fillet based on a Tth power 
polynomial and an exponential shape which closely resembled the above 
shapes. Inasmuch as no data are shown for the last two shapes the details 
of the fillet are restricted to the 9th power polynomial case. 

The shape of the fillet tangent to the flat plate at A and to the 
ramp at B where A and B are equidistant from C by the distance D is given 
by the 9th power polynomial 


y = 



The boundary conditions are: 
At point A 


y = 0 


2 3 4 

dy _ d^y _ d^y _ d^ 

dx 2 3 

dx dx dx 


H 


= 0 


(A-l) 


(A-2) 


At point B 
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y = D sin a 

^ = tan a 
dx 

..2 ,3 -.h 

= ±JL ^ ^ = Q 
2 3 4 

dx dx dx 

With the aid of the MIT algebraic computer language J4ACSYMA (ref. 21 ) the 
coefficients were found to be: 


(A-3) 
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where a is the ramp angle and D is the distance from the comer to 
fillet tangency points. The radius of curvature r is computed as 



(A-5) 


The profile of the fillet for D = .02 and a = 30° is shown as a solid 
line in figure 10. 
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APPENDIX B 


Incorporation of Velocity Profile into 
Integral Equation 


The boundary-layer integral equation vas derived in the main body of 
the paper as 


-0/_ n__ 

u 3u, /9u, u\f3u, 

\ - (ti: rrirj 

•'n ' Jr\ Jo 
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3 s 


r + n 2^ 
r R 


a /6 23 l \ + ^2 9u _ ^3^ 
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' (r + n 


n 
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-2 
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(r + n) 


dr , . 

^ — dri + V 
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^ — )>f(n) dn = 0 
van r+n/l 


(23) 


The assimed velocity profile that must be incorporated into the above 
integral equation is 


+ 

u = u^ [2.5 «<n (1 + n'^) + 5 .I - (3-387 n'^ + 5-l)e~°’ ] 


+ 0.5 Ug 




1 - cos 



(25) 
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In order to develop an expression for -r — note that u is a function of 

a s 

several variables that are in turn fimctions of s, hence, let 


u = [ug(s), u^(s), 6(s)] 


(B-1) 


or 



3u . 

9 u„ 9 s 

P 




9 u 9 6 
9 6 9 s 


Also, let 


or 


% ~ % 6(s)] 
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9 s 


9 u 9 u- 
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9 u 9 u 


9 u„ 9 s 
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Boundary conditions at n = 5 


u = 


u 

e 


(B-2) 


(B-3) 


(B-4) 


and 


9 u 

— ^=1 (B-5) 

^“6 


Substitution of equation (B-5) into equation (B-4), followed by the 

9ug 

elimination of from equation (B-2) with the use of the resultant 
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equation (B-h) yields 
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The derivatives 


9 u 


and 


96 


are not fimctions of n and it is convenient 


to remove them from the equation by differentiating equation (25) to give 
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(B-T) 


= u^|u^|r^ ~ ■ l-253196'^)e'°'^'^'^ j (B-8) 


so that equation (b- 6) can be written in the form 



9u 9 u \ ^ 

9 u ~ ^1 9 u_/ 9 s 
T 6 


/ 9 u 9 u \ 9 6 9 u ^ ^e 
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(B-9): 


9 U g g 9 u 

The remaining derivatives in the coefficients for -r , -r — , and -r in 

9 s 9 s 9 s 

equation (B-9) are functions of n. Expressions for these derivatives can 
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readily be deduced from equation (25) but are not reproduced here. 

In the region of the ramp fillet centrifugal forces will affect 

the boundary-layer edge velocity u^ if the radius of curvature exceeds 

the boundary-layer thickness 6. Consequently, the necessity that the 

edge velocities computed by the bomdary-layer equations match those 

computed from Prandtl-Meyer invicid flow requires that the unknown 
9 u 9 p 

6 w 

derivative -r be replaced by a new unknown -r . Put n = o in 

d S o S 

equations (21 ), (1T)j and (20) to obtain, respectively 
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(where the assumption was made that viscous effects at the edge of the 
boundary layer are small) 
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Elimination of and v from equation (B-ll) by use of equations (B-IO) 

os e 

and (3-11) results in 
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Substitution of the expression for from equation (B-9) into equation 

(B-13) yields after some simplification 
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Next, substitution of — from equation (B-2^) into equation (B-9) 
produces 
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Finally, the expression for -r — from equation (B-l6) is substituted into 

o S 

equations (23) and (2i+) for the strong interaction case and equation (23) 
for the weak interaction case to yield, after collection of like unknown 
terms, the algebraic equations 
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for the strong interaction and 
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for the weak interaction. 


The coefficients or functions E and Q are given in Appendix C. 

IJ J 


Ul 



APPENDIX C 
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The Functions E. . and Q. 
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The E..'s and Q.'s in equations (32) are as follows 
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and where , P^. P_, P. , P^ and , a^ have been defined in Appendix B. 
1 2 j 4 5 L c. 

these functions reduce to those of Tai (ref. 6). 
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:ure 1.- Theoretical model. 



Iiiviscid Solution for '.Vodge 
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Figure 2.- Effect of ramp angle. D = 0.02 
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Figure 3.- Comparison of calculated pressure distribution with experimental 
results. a = 10°, D = 0.02, (present calculations = 10X10®, 

R (experiment) = 2.1X106. 
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Figure 4.- Typical velocity profiles. D = 0.02, 
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Figure 5.- Effect of fillet length. a = 25° 
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Shang and Hankey (ref. 15) 
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